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ABSTRACT 

We study the brightness variations of galactic red supergiant stars using long-term visual 
light curves collected by the American Association of Variable Star Observers (AAVSO) over 
the last century. The full sample contains 48 red semiregular or irregular variable stars, with a 
mean time-span of observations of 61 years. We determine periods and period variability from 
analyses of power density spectra and time-frequency distributions. We find two significant 
periods in 18 stars. Most of these periods fall into two distinct groups, ranging from a few 
hundred to a few thousand days. Theoretical models imply fundamental, first and possibly 
second overtone mode pulsations for the shorter periods. Periods greater than 1000 days form 
a parallel period-luminosity relation that is similar to the Long Secondary Periods of the 
Asymptotic Giant Branch stars. A number of individual power spectra shows a single mode 
resolved into multiple peaks under a Lorentzian envelope, which we interpret as evidence for 
stochastic oscillations, presumably caused by the interplay of convection and pulsations. We 
find a strong 1 // noise component in the power spectra that is remarkably similar in almost all 
stars of the sample. This behaviour fits the picture of irregular photometric variability caused 
by large convection cells, analogous to the granulation background seen in the Sun. 
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1 INTRODUCTION 

Red supergiants (RSGs) are evolved, moderately massive (10-30 
M0) He-burning stars. As such, they are key agents of nucleosyn- 
thesis and chemical evolution of the Galaxy, and have also long 
been known for their slow optical variations. This variability is usu- 
ally attributed to radial pulsations (Stothers 1969, Wood, Bessell & 
Fox 1983, Heger et al. 1997, Guo & Li 2002), although irregu- 
lar variability caused by huge convection cells was also suggested 
from theory (Schwarzschild 1975, Antia et al. 1984) and observa- 
tions (e.g. Tuthill, Haniff & Baldwin 1997). The presence of os- 
cillations, in principle, offers the possibility of asteroseismology, 
i.e. measuring frequencies of stars and then identifying the modes 
of pulsations via comparing the observations to theoretical models. 
However, the time scales involved in RSGs makes this approach 
very difficult. Model calculations predict pulsational instability for 
fundamental and low-order overtone modes (Stothers 1972, Wood, 
Bessell & Fox 1983, Heger et al. 1997, Guo & Li 2002), with fun- 
damental periods ranging from 150 to 4,000 days. Consequently, 
any meaningful period determination may require decades of ob- 
servations. Here we analyse a homogeneous sample of RSG light 
curves with a typical time-span of over 22,000 days, which is the 
longest available observational data for these objects. 



Significant interest in RSG pulsations was driven by the dis- 
covery of their distinct period-luminosity (P-L) relation. Although 
the relation is not as tight as for Cepheids or red giant stars, vari- 
ous authors pointed out that because of their great intrinsic bright- 
nesses, RSGs might be useful as extragalactic distance indicators 
(e.g. Glass 1979, Feast et al. 1980, Wood & Bessell 1985). A re- 
cent example is that of Pierce et al. (2000), who re-calibrated the 
RSG P-L relation in the near-infrared and then measured the dis- 
tance to MlOl from 42 RSGs in that galaxy (Jurcevic et al. 2000). 
In this respect, galactic RSGs that are members of young 0-B as- 
sociations are particularly important because they can serve as local 
calibrators of the P-L relation - provided that reliable periods can 
be measured. Establishing this zero point is one of the by-products 
of this paper. 

Until recently, there was a serious discrepancy between evolu- 
tionary models and the empirical Hertzsprung-Russell diagram of 
RSGs. In particular, stellar evolutionary models seemed unable to 
produce RSGs that are as cool or as luminous as observed (Massey 
2003, Massey & Olsen 2003). A solution was found by Levesque et 
al. (2005), who derived a new (hotter) temperature scale from op- 
tical spectrophotometry for 74 Galactic RSGs, so that most of the 
discrepancy was removed. In addition, Massey et al. (2005) argued 
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that circumstellar dust around RSGs may account for many magni- 
tudes of extra extinction compared to other stars in the same 0-B 
associations, so that determining fundamental physical parameters 
can be quite difficult even for the brightest RSGs. The problem is 
very well illustrated, for instance, by the peculiar RSG VY CMa, 
for which previous radius estimates of up to 2,800 R© (Smith et 
al. 2001) were scaled down to 600 R© (Massey, Levesque & Plez 
2006), due to the hotter temperature scale alone. These wildly dif- 
ferent parameter values could be tested if pulsation mode identi- 
fication was made possible, which is another motivation for this 
work. 

In the General Catalogue of Variable Stars (GCVS, Kholopov 
et al. 1985-1988) there are two categories of variable RSGs: SRc 
and Lc. SRc type stars are semiregular late-type supergiants with 
amplitudes of about 1 mag and periods from 30 days to several 
thousand days. In contrast, Lc type stars are irregular variable su- 
pergiants having visual amplitudes of about 1 mag. As the distinc- 
tion between semiregular and irregular behaviour is not well de- 
fined in the literature, stars in both categories can possibly reveal 
important information. So far, the most detailed study of bright 
galactic RSG variables as a group was that by Stothers & Le- 
ung (1971), who discussed periods, luminosities and masses for 
22 stars. Their periods were either taken from papers published in 
the 1950s or (for 8 stars) determined by the authors using visual 
data. For some stars, such as a Ori, a Her and /i Cep, one can find 
more recent attempts to derive period(s) using various sources of 
photometric or radial velocity data (e.g. Mantegazza 1982, Smith, 
Patten & Goldberg 1989, Percy et al. 1996, Brelstaff et al. 1997, 
Rinehart et al. 2000). However, the typical time-span of analysed 
observations rarely exceeded 10-15 years and usually with many 
gaps, thus the derived periods (or mean cycle lengths) did not al- 
ways agree. 

The American Association of Variable Star Observers 
(AAVSO) has been collecting visual observations of variable stars 
for almost a century. The AAVSO database is by far the biggest 
one of its kind, containing over 12 million individual brightness es- 
timates for about 6,000 variable stars. Of these, roughly 50 stars 
are reasonably well-observed bright galactic RSGs of the GCVS 
types SRc and Lc. Using their observational records in the AAVSO 
database, we have carried out a period and light curve analysis of 
this group of stars. The main aim was to derive periods for all vari- 
ables. However, it turned out that in some cases the light curves had 
long enough time- span to study the long-term behaviour of their 
brightness fluctuations. Therefore, in addition to determining the 
dominant time-scales, we also discuss the nature and extent of the 
irregularities. 

The paper is organised as follows. In Sect. 2 we describe the 
sample selection, data processing and analysis. Sect. 3 contains the 
results, the newly determined periods; Sect. 4 discusses the multi- 
periodic nature and a comparison with pulsation model predictions. 
In Sect. 5 we present evidence for stochastic oscillations in some 
of the stars and the presence of a strong 1/f noise. Sect. 6 briefly 
summarizes the main findings of this study. 



2 SAMPLE SELECTION AND DATA ANALYSIS 

To find all well-observed RSG variables, we browsed through the 
AAVSO light curves of all stars listed as SRc or Lc type red su- 
pergiant in the GCVS. As an independent source of well-known 
galactic RSG stars, we checked the objects in Stothers & Leung 
(1972), Pierce et al. (2000) and Levesque et al. (2005). The latter 



paper was also used to collect the main physical parameters of the 
sample. The full set of stars is presented in Table [l] where, in ad- 
dition to the GCVS data, we also indicate the starting and ending 
Julian Dates of the AAVSO light curves and the number of points 
in the raw data. We kept stars with at least 400 individual points 
(actually, W Ind only has 399). In total, we retained 234,527 visual 
magnitude estimates for 48 stars, with a total time-span of 2,926 
years, or about 61 years of data for each star. 

The data were handled in a similar way to our previous anal- 
yses of long-term visual light curves (e.g. Kiss et al. 1999, 2000; 
Kiss & Szatmary 2002; Bedding et al. 2005). The raw hght curves 
were plotted and then inspected for outlying points, which were 
removed by a sigma-clipping procedure. 10-day bins were calcu- 
lated, which helped make the data distribution more even because 
in every star, there was gradual increase in the frequency of obser- 
vations by a factor of 3 to 5 over the 20th century. Without binning, 
any kind of period determination would have been strongly biased 
by the latter half of the data. 

Representative light curves are shown in Fig.[l] where the time 
axes have the same range to give a comparative illustration of the 
time- spans of the data. We see various characteristic features in the 
light curves. The most impressive stars are those with the largest 
amplitudes, like S Per and VX Sgr, where the full brightness range 
is similar to that of a Mira-type variable. However, none of the stars 
are as regular as a Mira. At the other end of the spectrum, we see 
stars with low amplitudes (<1 mag full range) and in many cases, 
two separate time- scales of variations: a slow one of a few thousand 
days and a faster one of a few hundred days (e.g. a Ori, TV Gem, 
a Her). In between, we also see intermediate- amplitude objects (of 
2 to 3 mags), whose light curves are dominated by relatively stable 
cycles, even suggesting the presence of stable oscillations (W Per, 
AH Sco). However, in many of the stars, there seems to be no stable 
periodicity, although the data clearly show real brightness fluctua- 
tions (e.g. /i Cep, BU Gem, CE Tau). 

To detect periodicities we calculated Fourier spectra of the 
light curves and then identified frequencies of power excesses. For 
this, we used Period04 of Lenz & Breger (2004) to carry out stan- 
dard iterative sine- wave fitting. In every iteration, a sine- wave cor- 
responding to the highest peak in the frequency spectrum was fit- 
ted and subtracted from the data. The spectrum was then recalcu- 
lated using the residual data. The iterative procedure was stopped 
when the residual spectrum did not contain peaks more than three 
times the noise floor. We then merged the frequencies in a few 
well-confined groups, within which the amplitude-weighted cen- 
troid defined a mean "cycle length" of the star. These were later 
refined using a quantile analysis (see below). In the stars with high- 
est amplitudes (e.g. S Per, VX Sgr) we found integer multiples of 
the dominant frequency, showing the non- sinusoidal shape of the 
curve; in those cases, the harmonics were not treated as separate 
periods. 

To allow direct comparison between the stars, and because 
the oscillation modes are resolved, we converted the power 
(=amplitude^) spectra to power density (power per frequency res- 
olution bin, expressed in mag^(c/d)~^), since the latter is indepen- 
dent of the length of the observations (Kjeldsen & Bedding 1995, 
App. A.l). Heavily smoothed power density spectra (PDS) were 
then used to determine the frequency dependence of the noise (ob- 
servational as well as astrophysical), which seem to be the domi- 
nant factor in some of the variables. 

Sample power spectra are plotted in Figs. 12131 where we also 
show the spectral window for each star. In high- amplitude stars 
close to the ecliptic plane (AH Sco, VX Sgr), we see strong yearly 
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Figure 1. Sample AAVSO light curves of red supergiant variables (10-day bins). Note the different magnitude scale in each plot. 
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Table 1. The basic properties of the sample. Classification and spectral types were taken from the General Catalog of Variable Stars; mmax and mmin are the 
maximum and minimum visual magnitudes as measured from the binned light curves. Every dataset ends in late 2005. 
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aliases, and that is why there are relatively high peaks in the spec- 
tra that remained unmarked. Also, for a few stars we find rela- 
tively strong peaks at a period of exactly 1-yr, whose reality is quite 
doubtful. Since these variables have very red colours, visual obser- 
vations can be affected by seasonal poor visibility, when the stars 
are observed at such high air-masses that the colour difference be- 
tween the variable and comparison stars may lead to differential 
extinction of a few tenths of a magnitude. We found very similar 
1-yr periods in semiregular red giants (Kiss et al. 1999), so that pe- 
riods between 355 and 375 days that were based on a single sharp 
peak were omitted from further analysis. 



As can be seen from the spectral windows in Figs. I2I3I the 
typical sampling is excellent and most of the structures in the PDS 
are real. The closely spaced peaks and their power distribution is 
very similar in shape to those of the solar-like oscillators and pul- 
sating red giants with stochastic behaviour (Bedding 2003, Bedding 
et al. 2005). For that reason, hereafter we make an important dis- 
tinction between the instantaneous period of the star, which can be 
measured from a shorter subset of the light curve and the mean pe- 
riod, which has the real physical meaning. Due to the seemingly 
irregular nature of the light variations, there are certain limitations 
in assigning "periods" to the observations. Assuming that stochas- 
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Figure 2. Power density spectra of AAVSO data. The insets show the spectral window on the same scale, while the arrows mark the adopted frequencies. 
Integer multiples of a main frequency (2/, 3/, etc.), yearly aliases (/+l/yr) or exactly 1-yr periods (+ signs) were not treated as separate periods. 



tic excitation and damping occur in these stars, one will always 
measure different period values from datasets that are comparable 
in length to the mode lifetime; however, that does not mean that 
the physical frequency of that particular mode has changed. In that 
sense we avoid using the term "period change", because what we 
can measure does not imply any change at all in the period of the 
pulsation mode. 

Further support to this was given by time-frequency analy- 
sis, which can reveal time-dependent modulations of the frequency 
content (e.g. Szatmary, Vinko & Gal 1994, Foster 1996, Bed- 



ding et al. 1998, Szatmary, Kiss & Bebesi 2003, Templeton, Mat- 
tel & Willson 2005). We checked several time-frequency distri- 
butions for gradual period evolution using the software package 
TIFRAN (Time FRequency ANalysis), developed by Z. Kollath 
and Z. Csubry at the Konkoly Observatory, Budapest (Kollath & 
Csubry 2006). In all cases, the wavelet maps clearly showed that 
close multiplets in the power spectra corresponded to a single peak 
that was subject to a time-dependent "jitter" with no signs of long- 
term evolution. Sample wavelet maps are shown in Fig. |4l with 
mean frequencies marked by the horizontal ticks. It is obvious, for 
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Figure 3. Sample power density spectra (cont.) 



example, that what seems to be two closely spaced peaks in AH Sco 
does in fact reflect a slight shift of of the dominant peak in the mid- 
dle of the dataset. Similarly, the broad power excess in W Per (Fig. 
13 corresponds to a peak that is highly variable in time. This be- 
haviour is very typical for most of the stars and is similar to that 
of the semiregular red giant variables (e.g. Percy et al 1996, 2003). 
The fluctuations of the instantaneous frequency usually do not ex- 
ceed 5-10% of the mean, although in extreme cases (like a Ori, 
TV Gem) the full width of a power excess hump can have relative 
width of about 20% in frequency. 

Keeping in mind these time-dependent changes, we deter- 



mined mean periods as follows. We identified probable dominant 
frequencies from the well-defined humps of power excesses in the 
PDS, as described above. Then we selected the low- and high- 
frequency limits of each hump as the frequencies where the power 
density reached the level of the noise. In this interval the cumula- 
tive distribution of the power was calculated and normalized to 1. 
The 0.5-quantile, the frequency where the normalized distribution 
crossed the 0.5 value is adopted as the mean frequency of the hump. 
Similarly, the width of the power excess was assigned to be the dif- 
ference between 0.166 and 0.834 quantiles, divided by 2. This way 
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Figure 4. Sample wavelet maps. The horizontal ticks on the right hand side 
show the mean frequencies. 



one can assume that at any given time the instantaneous period is 
within the [mean± width] interval with 1-sigma confidence. 

This method is somewhat subjective with the selection of the 
initial frequency intervals. But the quantiles are the most sensitive 
to the position of the peak and this stabilizes the results against 
the initial conditions. We tested this with comparing the results of 
the above discussed initial setup and another, strongly different test 
setup, where one of the initial positions was set twice further from 
the other one. The mean frequency was changed only by a few per 
cent, despite the obviously wrong initial settings. The width is al- 
most similarly stable with a variation of 5-20%. This showed that 
the uncertainty of the mean frequency is usually in the order of 
a few percent, being much smaller than the natural jitter present 
in the stars. In several cases we fitted Lorentzian envelopes to the 
power excess humps (see details in Sect. 5) to measure the centroid 
frequency and its damping rate. For every star with well-defined, 
regular humps the results were very similar to those of the quantile 
analysis. However, for other stars the PDS is quite noisy and the as- 
sumption of Lorentzian power distribution is not true, so that fitting 
Lorentzians was not possible for the whole sample. 
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Figure 5. A comparison of 33 periods for 28 stars. 



3 RESULTS 

We present the derived periods in Table |2| In total, we determined 
56 periods for 37 stars; in 6 cases ours is the first period determina- 
tion in the literature. We found published periods for 31 stars, but 
except for a few well-studies variables (e.g. /x Cep, a Ori, S Per), 
most of them were quite neglected in the past three decades. Never- 
theless, the overall agreement with the catalogued periods is good: 
for the shortest period stars (e.g. SS And, T Cet, Y Lyn, W Tri) our 
findings are in perfect agreement with the GCVS or other, more 
recent studies (with the differences staying below 2-3%). In a few 
cases we could not infer any clear coherent signal from the curves, 
so we show the upper limits on the semi- amplitude in parentheses. 

In Fig.|5]we compare our periods with the literature. For this 
comparison we selected those periods from the literature that were 
definitely corresponding to ours. There are three significant out- 
liers, namely CL Car at 2600 d (with 952 d from ASAS project), 
ST Cep at 3300 d (with 2050 d from Stothers & Leung 1971) and 
XX Per at 3100 d (with 4100 d also from Stothers & Leung 1971). 
In case of CL Car, we checked the ASAS observations of the star 
and found that the 952 d period is an error. For ST Cep and XX Per 
our data have significantly longer time- span, so that the newly de- 
termined values are likely to be more accurate. 

We have attempted to confirm the red supergiant status of 
stars in our sample. The minimum absolute bolometric magnitude 
for red supergiants is about —5 mag (Meynet & Maeder 2003), 
which can be translated to Mk ~ —8 mag using the relation 
mboi ~ + 3 that connects extinction-corrected bolometric and 
iT-band magnitudes (Josselin et al. 2000). In our sample, bolomet- 
ric magnitudes were recently determined for 18 stars by Levesque 
et al. (2005) and of these only ST Cep is located near the minimum 
luminosity, with Mboi = —5.48 mag. A few stars have useful Hip- 
parcos parallaxes {a Ori, a Sco, a Her, CE Tau), and all are at least 
1 mag brighter than Mk = —8 mag. A particularly interesting 
star is a Her because its period (125 d) is very short for a super- 
giant star. For instance. Jura & Kleinmann (1990) rejected all red 
supergiant candidates with periods less than 150 d, although it was 
clearly established from long-period variables in the Large Magel- 
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Table 2. Periods from this study and the Hterature. The given ±range values were calculated from the width of the power distribution in the frequency spectra 
and are dominated by the intrinsic jitter of the stars. Numbers in parentheses refer to an upper limit of amplitude (in magnitudes) in the Fourier spectrum in 
those cases where there was no peak with S/N>3, while (1/f) refers to cases where the spectra show a constant rise towards the lowest frequencies with no 
discernible peak. Stars in parentheses were rejected as supergiants. 
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lanic Cloud that the least luminous red supergiants can have periods 
between 100 and 150 d (Wood, Bessell & Fox 1983). For a Her, 
K = -3.70 mag (Richichi & Percheron 2002) and tt = 8.5 ± 2.8 
mas (ESA 1997), implying Mk = — 9.0 =b 0.7 mag, so our con- 
clusion is that shorter periods (100 — 150 d) can indeed exist in 
RSGs. 

For four stars, whose names are shown in parentheses in Ta- 
ble |2l the RSG class is somewhat doubtful. For UZ CMa, both the 
AAVSO light curve and the ASAS F-band observations infer a pe- 
riodicity around 40 days; the ASAS CCD-V light curve is very 
typical of a short period red giant semiregular variable, presumably 
on the AGB, which is also supported by the luminosity class 11. 
IS Gem, although classified as SRc, has too early a spectral type 
(K3II) and its i^-band absolute magnitude is only —2.10 mag, be- 
ing consistent with the luminosity class. For T Cet, the periods and 



their ratio of 1.87 are fairly typical for semiregular AGB stars (Kiss 
et al. 1999), and its Hipparcos parallax and 2MASS K magnitude 
imply Mk = — 7.7 ± 0.4 mag, which is too faint. Y Lyn is a very 
typical semiregular variable with a long- secondary period, which is 
a well-known but still poorly understood phenomenon in AGB stars 
(Wood, Olivier & Kawaler 2004). Moreover, its i^-band absolute 
magnitude from the Hipparcos parallax and 2MASS i^-magnitude 
is only —7.7 ± 0.7 mag. For these reasons, hereafter we exclude 
these four stars from further analysis. 



4 DISCUSSION OF THE MULTIPERIODIC NATURE 

According to the theoretical predictions, multiple periodicity may 
arise from multimode pulsations. Radial oscillations have long 
been predicted by model calculations, starting from the pioneering 
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Figure 6. Period-absolute magnitude relations for 30 periods of 18 stars. 
Triangles refer to RSG variables in the LMC, with periods and photometry 
taken from Wood, Bessell & Fox (1983). Solid circles and asterisks were 
divided at P=1000 d. The dotted lines are based on model calculations by 
Guo & Li (2002) for solar metallicity, while the solid lines shows the best-fit 
period-K magnitude relation. 

work of Stothers (1969). Two extensive investigations on pulsation 
properties of RSGs were recently published by Heger et al. (1997) 
and Guo & Li (2002). Although they treated convection in a differ- 
ent way (Heger et al. adopted the Ledoux criterion and treated semi- 
convection according to Langer, Sugimoto & Fricke (1983); Guo 
& Li used the mixing-length theory of Bohm-Vitense (1958) and 
adopted the Schwarzschild criterion to determine the boundaries of 
convection and semiconvection zones), both studies confirmed ear- 
lier results and predicted excitation of the fundamental, first and 
possibly second overtone modes. The fundamental mode's growth 
rate always exceeded that of the overtone modes and was found 
to increase with luminosity for a given mass. On the other hand, 
studies of convection in the envelopes of red giants and supergiants 
showed that the dominant convective elements can be comparable 
in size to the stellar radius, which could explain both the observed 
irregular variations (Schwarzschild 1975, Antia et al. 1984) and 
the hotspots on late- type supergiants, detected using interferomet- 
ric techniques (Tuthill, Haniff, & Baldwin 1997). A similar mech- 
anism was proposed by Stothers & Leung (1971) for explaining 
the longer periods of RSG variables, arguing that the similarity of 
those periods and the time scale of convective turnover suggests a 
link between the two phenomena. Interestingly, periods that are too 
long were also found in other luminous stars. Maeder (1980) cal- 
culated empirical pulsation constants of blue-yellow supergiants, 
which were found to be systematically larger than the theoretical 
Q value for the fundamental mode of radial oscillation. Maeder 
(1980) suggested that those long periods can be due to non-radial 
oscillations of gravity modes. The same mechanism was proposed 
by Wood, Olivier & Kawaler (2004) as a possible explanation of 
the long- secondary periods of AGB stars. 

4.1 Period-luminosity relations 

A significant fraction of our sample can be characterized by two 
dominant periods, one of a few hundred days and one of 1500-2000 
days. The strong similarity that we found in many cases argues for 
the reality of these long periods. Some of the light curves are more 
than twice as long as those analysed by Stother & Leung (1971) 



and the fact that we derive very similar values for the long peri- 
ods indicates the reliability of the results. To reveal deeper insights 
into the nature of the multiply periodic variations, we studied the 
period-luminosity distribution for all stars having a useful estimate 
of luminosity. 

To construct the period-iC-band absolute magnitude diagram, 
we took the following steps: {i) Mboi values from Levesque et al. 
(2005) were converted to Mk using the Josselin et al. (2000) re- 
lation; { ii) for a Her, a Ori and CE Tau we used Hipparcos paral- 
laxes and K-band data to calculate Mk directly; (Hi) for compar- 
ison we added RSGs in the Large Magellanic Cloud, taken from 
Wood, Bessell & Fox (1983). The resulting P-L diagram is shown 
in Fig.|6l where the lines represent fundamental, first overtone and 
second overtone modes for models of solar metallicity by Guo & 
Li (2002). The help distinguish between the two sequences, periods 
longer than 1000 d were plotted with different symbols. 

The shorter periods are well matched by the fundamental and 
first overtone modes of the models. Moreover, there is perfect 
agreement with the LMC red supergiant sample, too. /a Cep and 
PZ Cas are above the luminosity range of the models and it is pos- 
sible that in those cases the longer period is the fundamental mode. 
RT Car and maybe IX Car seem to be too luminous for fundamen- 
tal pulsation, but otherwise models strongly favour fundamental or 
first overtone modes for the shorter periods. It is worth noting that 
a Her has an almost identical counterpart in the LMC, which pro- 
vides a retrospective confirmation of the reality of the 125 d period. 

Excluding the short period of RT Car, we fitted the following 
period-iT magnitude relation to the solid dots in Fig.|6l 

Mk = (-3.44 ±0.6) log P + (-1.6 ± 1.6) 
with 0.46 mag rms, which is about the typical uncertainty of indi- 
vidual absolute magnitudes. The slope of this relation is very sim- 
ilar to the slope of the Mira P-L relation in the Large Magellanic 
Cloud (-3.52 ± 0.03, Ita et al. 2004). The agreement shows the 
similarity of pulsations in red giant and red supergiant stars. How- 
ever, the zero points are very different, as the Ita et al. relation for 
fundamental mode AGB stars has a zero point of -i-1.54±0.08 mag 
(with fiLMC — 18.50). 

The long secondary periods (LSPs) remain beyond the limits 
of the models, and both the updated models and the refined phys- 
ical parameters of the stars keep the original conclusion on the 
peculiar nature of LSPs by Stothers & Leung (1971) unchanged. 
They cannot be explained by radial pulsations, because the period 
of the fundamental mode is the longest possible for that kind of 
oscillation. Metallicity effects also cannot explain the LSPs. Look- 
ing at the positions, for instance, a Her or a Ori, their absolute 
magnitudes are 1-2 mags fainter than those of the solar metallic- 
ity models. According to Guo & Li (2002), ^Mboi ^ 0.835 log Z, 
which means that unphysically large metallicities (10-100 times 
solar) would be needed to account for the low luminosities. A pos- 
sibility is that heavy circumstellar extinction in K-bmid makes the 
stars fainter. Massey et al. (2005) indeed found many magnitudes 
of circumstellar F-band extinction in a number of galactic RSGs, 
which showed that the effect may not be negligible (for example, 
in Fig. |6la Her, a Ori and CE Tau were not corrected for this). 
Furthermore, it is known that the reddest Mira stars in the Small 
Magellanic Cloud have K-bmid magnitudes that are fainter by 1- 
2 mag than predicted by the Mira P-L relation (Kiss & Bedding 
2004), which is similar to what we see here for the red supergiants. 
On the other hand, the very good agreement for the shorter pe- 
riods between the galactic and the LMC samples argues against 
this explanation (see, e.g., a Her and the closest LMC point in 
Fig.|6j. Also, the 4-5 mags extra optical extinction that has been 
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Figure 7. Period-X -magnitude relations for the galactic sample (black cir- 
cles), RSG variables in the LMC (red triangles) and red giant variables from 
the MACHO database (blue dots, data taken from Derekas et al. 2006, as- 
suming )LiLMC = 18.50). Labels A, B, C and D were adopted from Wood 
(2000). 



found in cluster RSGs by Massey et al. (2005) could hardly explain 
the 1-2 mags extra extinction in the ii^-band, unless the extinc- 
tion law of the circumstellar matter is extremely different of the 
"standard" one {Ak/Av = 0.112, Schlegel, Finkbeiner & Davis 
1998). Therefore, the LSP P-L sequence is a separate entity, whose 
origins needs further investigation (possibly together with the AGB 
LSP phenomenon). 

In Fig. we compare the period-i^-magnitude relations for 
RSGs in our galactic sample (large black dots) and in the LMC (red 
triangles) with red- giant variables in the LMC (small blue dots), 
as observed by the MACHO project (data taken from Derekas et 
al. 2006). The MACHO sample comprises low-mass stars on the 
RGB and AGB. There is some similarity between the P-L rela- 
tions of the supergiants and the less luminous RGB and AGB stars, 
although the lack of precise distance estimates for our sample of 
RSGs makes a detailed comparison difficult. Also, we expect the 
RSG sequences to be broad because these stars cover a large range 
of stellar masses. Nevertheless, it does appear that the shorter pe- 
riods of the RSG stars mostly align better with Sequence B of the 
low-mass stars rather than with Sequence C, which would imply 
pulsation in the first-overtone rather than the fundamental. It would 
clearly be valuable to have long time series for a large sample of 
LMC stars, which could be obtained by analysing the photographic 
plate archives. 



4.2 Pulsation constants 

Further interesting details are revealed by examining the pul- 
sation constants Q = P{M /MqY^^ {R/ Rq)-^^^ and W = 
P{M/Mq)(R/Rq)~'^. The first is the classical period-density re- 
lation, while the second is the natural form of the pulsation con- 
stant if the oscillations are confined to the upper layers of the en- 
velope (Gough, Ostriker & Stobie 1965). Stothers (1972) showed 
that for pulsation models of massive red supergiants. Wo of the 
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Figure 8. Pulsation constant Q and W. The dashed and dotted lines show 
the fundamental and first overtone models of Guo & Li (2002). The down- 
ward outlier at P = 200 d is RT Car. 



fundamental mode was more constant than Qo. We have calcu- 
lated both quantities for the stars in Fig. |6l as follows. Radius 
values were either taken from Levesque et al. (2005) or calcu- 
lated from the parallax-based luminosity (a Her and CE Tau) and 
temperature (a Her: Levesque et al. 2005; CE Tau: Wasatonic & 
Guinan 1998). Masses were estimated from the approximate rela- 
tion log(M/M0) = 0.50 - O.lOMboi, which comes from evolu- 
tionary calculations (Levesque et al. 2005); for this, Mboi values 
were taken from Levesque et al. (2005) or calculated from Mk. 

We plot the resulting pulsation constants Q and W as func- 
tion of period in Fig. |8] In addition to the empirical values, we 
also show fundamental- and first overtone-mode models of Guo 
& Li (2002). Based on this diagram we draw several conclusions. 
Firstly, looking at the model calculations, W of the fundamen- 
tal mode is indeed a better constant than Q, in accordance with 
the theoretical expectations, while the first overtone-mode models 
show an opposite behaviour. Secondly, all periods less than 1000 
d agree with fundamental pulsation for both Q and W (except 
for RT Car). Thirdly, there is a remarkable feature of the pulsa- 
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tion constants: the mean value of Q for the shorter periods and 
the mean value of W for the long periods are practically equal: 
(Q) = 0.085 ± 0.02 and {W) = 0.082 zb 0.03. Although this 
might be coincidence, it may suggest that the fundamental period 
and the long secondary period are intimately connected via the rela- 
tion Po/Plsp = {M/MqY^'^{R/Rq)-^^^, which could be used 
as test for future models of the LSPs. 

One application of the two pulsation constants could be check- 
ing the consistency of a given set of physical parameters for a 
star with the observed periodicity. A notoriously ambiguous case is 
VY CMa, for which vastly different radius estimates can be found 
in the literature, ranging from 600 Rq up to 2,800 R© (Massey, 
Levesque & Plez 2006 and references therein). For example, Mon- 
nier et al. (1999) adopted M ^ 25 M© and R ^ 2,000 R©, 
while Massey, Levesque & Plez (2006), using the new temper- 
ature scale of Levesque et al. (2005), arrived to M ^ 15 M© 
and R ^ 600 R©. Using the mean cycle length of 1,600 d, the 
high-mass/large radius parameter set results in Q = 0.089 and 
W — 0.01, both being consistent with fundamental mode pul- 
sation. For the low-mass/small radius set, the two constants are 
Q = 0.42 and W = 0.066, which could be acceptable if the 1600 
d period referred to a long secondary period. At this stage, unfor- 
tunately, the periodicity does not help solve the problem, but since 
no other RSG has a fundamental-mode period greater than 1000 d, 
we have a slight preference for the low-mass/small radius set and 
the LSP interpretation of the 1,600 d periodicity. 



5 EVIDENCE FOR STOCHASTIC OSCILLATION AND 
l/F NOISE 

5.1 Lorentzian envelopes in the power spectra 

The AAVSO observations span many decades and sometimes al- 
most a century, thus we are now in a better position to understand 
the irregularity of RSG variability than previous researchers. Not 
much effort was put into this direction in the past: irregularity was 
taken as a general description of non-periodic brightness fluctua- 
tion in red giants and supergiants. For semiregular red giants the 
nature of irregularities was addressed by Lebzelter, Kiss & Hinkle 
(2000), who compared simultaneous light and velocity variations 
in a sample of late-type semiregular variables (of the GCVS types 
SRa and SRb) and concluded that the observed variability is most 
likely a combination of pulsations and additional irregularity in- 
troduced by, e.g., large convective cells. The latter phenomenon is 
theoretically expected (Schwarzschild 1975, Antia et al. 1984), al- 
though not much is known on the time-dependent behaviour of the 
integrated flux variations that arise from the huge convective cells. 
Percy et al. (2003) found growth/decay timescales of 1-5 years in 
a sample of small- amplitude pulsating red giant stars, which they 
interpreted as the natural growth (or decay) times for the pulsation 
modes. Recently, Bedding et al. (2005) discussed the possibility 
of solar-like excitation of the semiregular variable L2 Pup. In that 
star the the power distribution closely resembles that of a stochasti- 
cally excited damped oscillator and is strikingly similar to close-up 
views of individual peaks in the power spectrum of solar oscilla- 
tions: there is a single mode in the spectrum which is resolved into 
multiple peaks under a well-defined Lorentzian envelope. The en- 
velope's width gives the damping time (or mode lifetime), which is 
one of the main characteristics of a damped oscillator. Bedding et 
al. (2005) argued that the close similarity may imply solar-like ex- 
citation of oscillations, presumably driven by convection. Another 
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Figure 9. Lorentzian fits (thick lines) of the power spectra (thin lines). 



interpretation could be that the cavity of oscillations is changing 
stochastically due to the convective motions, which affects the reg- 
ularity of the pulsations. 

We find the very same structures in the power spectra of sev- 
eral RSGs: a well-defined Lorentzian envelope under which the 
power is split into a series of narrow peaks. Probably the best ex- 
ample is W Per (Fig.|3}, where the window function is practically 
free of any kind of alias structure, so that no false structure due 
to poor sampling arises in the power spectrum. Other convincing 
examples include a Ori, TV Gem, S Per, AH Sco, VY CMa and 
a Her, whereas in a few stars the power is so spread over a large 
range of frequencies that no regular envelope can be traced in the 
spectrum. 

Adopting the same approach as Bedding et al. (2005), i.e. as- 
suming a stochastically excited damped oscillator, we have fitted 
Lorentzian profiles to the power spectra assuming statistics with 
two degrees of freedom. This is based on a maximum-likelihood fit, 
assuming an exponential distribution of the noise (Anderson, Du- 
vall & Jefferies 1990, Toutain & FroWich 1992). The fit gives the 
centroid frequency and the half-width at half power F, which can 
be converted to mode lifetime via r = (27rr)~^. 

We show the clearest examples of the Lorentzian fits in Fig.|9l 
We have looked for correlations between mode lifetime, luminos- 
ity and length of the long secondary periods. Generally, the mode 
lifetime is typically 3-4 times the pulsation cycle, with notable ex- 
ceptions of TV Gem (less than 1) and a Her (about 8). There might 
be a slight correlation between the long secondary period and the 
mode lifetime, but the small number of multiperiodic stars with 
well-defined Lorentzian fits prevented a firm conclusion. 
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5.2 1/f noise 

We also investigated the nature of irregularities through the shape 
of the noise level in the power spectra. The analysis of fluctua- 
tion power spectra density is a common tool in studies of unpre- 
dictable and seemingly aperiodic variability, i.e. noise that arises 
from a stochastic process. In this context, the noise is intrinsic to 
the source and not a result of measurement errors (such as Pois- 
son noise). In astrophysics. X-ray light curves of active galaxies 
and interacting binaries have been a major inspiration of such stud- 
ies (e.g. Vaughan et al. 2003 and references therein). Of particu- 
lar interest are noise series whose power spectra are inverse power 
functions of frequency, the so-called noises, for which exam- 
ples have been found in a wide range of natural phenomena (Press 
1978). Bak, Tang & Wiesenfeld (1987) showed that dynamical sys- 
tems with spatial degrees of freedom naturally evolve into a self- 
organized point. 1/f noise in these systems simply reflects the dy- 
namics of a self-organized critical state of minimally stable clusters 
of all length scales, which in turn generates fluctuations on all time 
scales. Specifically, turbulence is a phenomenon, for which self- 
similar scaling is believed to occur both in time and space. There- 
fore, combining the Bak, Tang & Wiesenfeld (1987) theory with the 
Schwarzschild (1975) mechanism of producing irregular variability 
via large convective cells, one expects a strong 1/f noise compo- 
nent in the power spectra of RSG brightness fluctuatios, similarly 
to solar granulation background (e.g. Rabello-Soares et al. 1997). 

That is exactly what we find for the vast majority of the sam- 
ple. In Fig. [To] we plot smoothed power density spectra in log-log 
representation. These were calculated from the initial spectra af- 
ter transforming to the log-log scale and then binning them with a 
stepsize of 0.1 dex. Besides 17 RSGs, we also show the averaged 
power density spectrum of the AAVSO observations of the classical 
Cepheid variable X Cygni, which varies between 5.9-6.9 
mag with a period of 16.38 d. This star, being a strictly regular pul- 



sating variable, served as a test object for the power distribution of 
visual observational errors. As indicated by the lower right panel 
of Fig.[T0l the power distribution is much flatter for almost three 
orders of magnitude, with only a slight rise towards the lowest fre- 
quencies. This is likely to arise from the fact that not many AAVSO 
observers are active for more than 10,000 days, so there is a fluctu- 
ation of observers on time-scales of longer than a few 1000 days. It 
is also apparent that the mean power in X Cyg is much lower than 
in any of the RSGs, scattering between log PD = 0-1. Its bright- 
ness range is roughly in the middle of the sample (see Table[I), the 
time- span is about 20,000 days, so that the mean noise level should 
be representative of the typical AAVSO light curve in our set. 

The RSG spectra look remarkably similar to each other. There 
is a roughly linear power increase towards the low-frequency end, 
which has very similar slope in almost all cases. Generally, we 
do not detect any flattening at the lowest frequencies, except for 
/ < 1/Tobs (Tohs is the time-span of the data), which suggests 
that the observations are too short to extend over the whole range of 
time- scales of fluctuations that is present in these stars (with a Her 
being an exception). Formal linear fits of the spectra resulted in 
slopes ranging from 0.8 to 1.2, but none was significantly different 
of 1 . Furthermore, distinct features in the spectra (like the peaks of 
the periods) made it difficult to fit the overall slope of the "noise 
continuum" accurately, and because of that we fixed the slopes at 
1 and fitted the zero-points only, excluding the frequency ranges of 
the higher power concentrations. The results are indicated by the 
thick black lines in Fig. 1101 

Based on the close similarity of the log-log power spectra we 
conclude that there is a universal frequency scaling behaviour in the 
brightness fluctuations of red supergiant stars that fits very well the 
expectations for background noise from convection. It is not sur- 
prising that period determination is so difficult for these variables: 
the longer we observe, the more power will be detected in the low- 
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frequency range, suggesting the presence of more and more "pe- 
riods" that are comparable to the full length of the data. Periods 
identified with the fundamental mode are definitely real - although 
strongly affected by the high-noise in the systems. We think most of 
the long- secondary periods are also real, because we see a consis- 
tent picture in many different stars. However, random peaks at the 
lowest frequencies naturally develop due to the strong 1// noise 
signal and it would be misleading to interpret each peaks as peri- 
ods. Particularly good examples are AD Per and W Cep, of which 
the latter has variations up to 2 mags in visual but the power fol- 
lows a well-defined 1//^ distribution (with a slope that may be a 
bit larger than 1). The whole set of phenomena is strikingly similar 
to what is observed in the Sun and other solar-like oscillators. The 
photometric granulation noise in a main- sequence solar-like oscil- 
lator has a time- scale of minutes and micromagnitude amplitudes. 
At the other end of the spectrum we see these pulsating red su- 
pergiants with time-scales of years and noise amplitudes reaching 
tenths of a magnitude; the underlying physical mechanisms seem 
to be the same all across the Hertzsprung-Russell-diagram. 

Finally, there is interesting correlation between the zeropoints 
of the fitted lines and the light curve amplitudes. The two highest 
amplitude stars, S Per and VX Sgr, have the highest noise levels, 
reaching one to two orders of magnitudes higher at log / < — 4 
than any other variable. This shows they are fundamentally more 
dynamic in pulsation and noise generation, which reminds us the 
predicted "superwind" phase of RSGs just preceding the supernova 
explosion. Pulsation model calculations by Heger et al. (1997) have 
shown that very large pulsation periods, amplitudes and mass-loss 
rates may be expected to occur at and beyond central helium ex- 
haustion over the time-scale of the last few 10^ years. The physical 
reason for this is the resonant character of pulsations when the pul- 
sation period and the Kelvin-Helmholtz time-scale of the pulsating 
envelope evolve into the same order of magnitude. A similar re- 
sult was found by Bono & Panagia (2000), whose pulsation models 
showed larger amplitude and more irregular variations for lower 
values of Teff, which turned out to be the main governer of the pul- 
sational behaviour. As Lekht et al. (2005) noted for S Per, overall 
dimming of the star after a period of stronger oscillations may be 
due to subsequent enhanced mass-loss and ejection of a dust shell 
that screens the stellar radiation. This is clearly seen for S Per and 
VX Sgr, and maybe in AH Sco and VY CMa, giving supporting 
evidence for variations of the pulsation driven mass-loss on a time- 
scale of ~20 years. 



6 CONCLUSIONS 

Red supergiants as variable stars have long been known for their 
semiregular brightness fluctuations. Using the extremely valuable 
database of visual observations of the AAVSO, we were able to 
study the main characteristics of their regular and irregular varia- 
tions. From a detailed analysis of power spectra and time-frequency 
distributions, this paper discussed the properties of pulsations and 
their physical implications. 

The sample contains several types of light variations. A few 
stars (S Per, VX Sgr, AH Sco) have very large peak-to-peak ampli- 
tudes that may reach up to 4 magnitudes in the visual. These objects 
show both the most coherent periodic signals and the highest level 
of 1// noise in the power spectrum, which might be a sign of a "su- 
perwind" phase of RSGs that precedes the supernova explosion. We 
argue that multiple periods found by other studies were sometimes 
artifacts caused by the extreme noise levels of the systems. A more 



common type of variability is characterized by two distinct peri- 
ods, one of a few hundred days and one of a few thousand days. 
The archetypes of these stars are a Ori and TV Gem, both having a 
period around 400 days and another around 2000 days. The shorter 
periods can be identified with the radial fundamental or low-order 
overtone modes of pulsation, while the longer one is very similar 
to the Long Secondary Periods of AGB stars, whose origin is not 
known yet, but could be due to binarity, magnetic activity or non- 
radial ^-modes (Wood et al. 2004). Besides fundamental pulsation, 
we also see evidence for first and possibly second overtone modes. 
Finally, in a few stars we cannot infer any periodicity at all (e.g. 
W Cep, AD Per) and these are the truly irregular variables domi- 
nated by the 1// noise in the power spectrum. 

The period jitter of the pulsation modes produces in the power 
spectrum a well-defined Lorentzian envelope. Interpreting this as 
evidence for stochastically excited and damped oscillations, we 
measured the mode-lifetime (or damping rate) in red supergiants 
for the first time. In most stars it is several times longer than the 
period of pulsations. Since the damping rate depends on the stellar 
structure and convection properties in a complex way, with many 
weakly constrained theoretical parameters (e.g. Balmforth 1992) 
and, moreover, currently there are no theoretical calculations di- 
rectly applicable to red supergiants, it is impossible to qualify the 
agreement with theoretical expectations. However, the strong 1// 
noise component that seems to be ubiquitous in the whole sample, 
strongly favours the Schwarzschild (1975) mechanism of produc- 
ing random brightness variations with huge convection cells, anal- 
ogous to the granulation background seen in the Sun. 

Finally, stochastic oscillations discussed in this paper may of- 
fer an explanation for the seemingly random behaviour of less lu- 
minous red giants stars. For example, random cycle-to-cycle fluc- 
tuations of the periods of Mira stars (e.g. Eddington & Plakidis 
1929, Percy & Colivas 1999) are likely to be analogous to the pe- 
riod jitter we found here. It would be helpful to have theoretical 
models that are specifically designed for these phenomena, because 
the extensive records of homogeneous visual observations allow us 
measuring noise properties of red giant and supergiant stars quite 
accurately. In principle, internal physics could be probed through 
these accurate measurements provided that realistic models are cal- 
culated. 
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